Spin-polarization coupling in multiferroic transition-metal oxides 
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A systematic microscopic theory of magnetically induced ferroelectricity and lattice modulation is 
presented for all electron configurations of Mott-insulating transition-metal oxides. Various mecha- 
nisms of polarization are identified in terms of a strong-coupling perturbation theory. Especially, the 
spin-orbit interaction acting on the ligand p orbitals is shown to give the ferroelectric polarization of 
the spin-current form, which plays a crucial role particularly in e g systems. Semiquantitative agree- 
ments with the multiferroic TbMnOs are obtained. Predictions for X-ray and neutron scattering 
experiments are proposed to clarify the microscopic mechanism of the spin-polarization coupling in 
different materials. 

PACS numbers: 75.80.+q, 71.70.Ej, 77.80.-e 



The coupling among the charge, spin, and orbital de- 
grees of freedom in Mott insulators has been one of the 
central issues in strongly correlated electron systems. De- 
spite the charge localization, its displacement can occur 
accompanied by the lattice deformation as in the dielec- 
tric or ferroelectric band insulators. In the case of Mott 
insulators, the charge displacement is usually triggered 
by a non-trivial low-energy spin and/or orbital structure, 
which is unveiled within an energy scale of the exchange 
interaction J ~ V^ 4 /A 3 with the hybridization V and the 
charge transfer energy A between the transition-metal 
(TM) d and ligand (L) p orbitals. There, the ferroelec- 
tricity and the magnetism sometimes appear concomi- 
tantly as a spontaneous magnetoelectric effect dubbed 
multiferroic order [1-3]. Usually, the largest magneto- 
electric force is given by the exchange-striction. Spatial 
imbalance in the symmetric magnetic exchange leads to 
a lattice modulation. When it breaks the inversion sym- 
metry, a finite net dipole moments appears as the ferro- 
electric polarization. 

Recent observations of the ferroelectric polarization at 
the onset of the spiral magnetic order in TbMnC>3 [4-6] 
and M3V2O8 [7] have attracted revived interests in multi- 
ferroic behavior from the viewpoints of the fundamental 
physics as well as applications [3]. Microscopic [8] and 
phenomenological [9, 10] theories have revealed that spi- 
ral and conical spin structures can couple to the polariza- 
tion P in the absence of any symmetry breaking through 
P J s , where J s oc ex(S r xS r+e ) oc S(V-S)-(S- V)S is 
the spin current produced between the two spins S r and 
S r+e at sites r and r + e, respectively. This coupling has 
been derived by taking account of the spin-orbit coupling 
A for a model of the ti g orbitals of TM ions hybridized 
with the L p orbitals [8, 11], an e g system coupled to the 
localized tig spin [12], and a sp hybridized model [13]. 



Effects of the exchange-striction and the Dzyaloshinskii- 
Moriya (DM) interaction have also been discussed in the 
context of the multiferroic i?MnC>3 [14]. 

In view of the recent proliferation in the number 
of multiferroic materials, a proper classification scheme 
seems to be in order. In this Letter, we propose a 
systematic microscopic theory of the magnetically in- 
duced polarization and lattice modulation in multiferroic 
transition-metal oxides, by means of a strong-coupling 
expansion in both V/A and A/ A. We classify various 
microscopic mechanisms for the polarization and their 
distinct features. The classification scheme offers a way 
of identifying the material-specific mechanism by prob- 
ing the magnetically induced non-uniform shifts of the L 
ions. As the specific demonstration, key experimental ob- 
servations on the polarization and the lattice modulation 
in TbMnC>3 are analyzed. 

Let us consider a cluster composed of TM ions at r and 
r + e hybridized through a L ion at the center r + e/2 as 
shown in Fig. 1 (a). For simplicity, we neglect a deviation 
of the bond angle away from 180°. The dipole moment 
Pr+e/2 induced at the ligand site r + e/2 is generally 
expressed as 



r> _ [ pnm 

*>+e/2 — l/r.r+e 



; (m T . 



3)] 



+P™r+e I( e ' rn r )m r - (e ■ m r+e )m r+e ] 
+Pr P r+e e x (m r x m r+e ) (1) 

up to second order in the local spin moments m r and 
m r+e . Here, P nm oc V/A is the nonmagnetic term that 
occurs regardless of the magnetic properties. P ms oc 
(V/A) 3 represents the magnetostriction term. These ef- 
fects do not demand the spin-orbit interaction. The or- 
bital term P orb ~ mm(A/F, l)(V/A) has been obtained 
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FIG. 1: (color online) (a) The TM-L-TM cluster model with 
the rod-type d 3x 2_ r 2/d 3y 2_ r 2 staggered orbital order under 
a noncollinear spin configuration with the associated electric 
polarization, (b) The lattice structure of the perovskite sys- 
tem within the xy (or ab) plane. The staggered orbital order 
for TbMnOa is also shown, (c) The level scheme for the i| g eg 
high-spin (S = 2) configuration. 



in the study of partially filled degenerate t2 9 orbitalsfll]. 
Finally, the spin-current term P sp oc (A/A)(V/A) 3 was 
demonstrated for the degenerate ti g model [8, 11] and, 
as will be shown below, exists for any Mott-insulating 
electron configuration. P sp and P orb are even under the 
spatial inversion and found only in the presence of spin- 
orbit interaction, while P nm and P ms are odd and present 
only in the absence of the inversion symmetry. 

On an array of TM-L-TM-L with (at 

most) two inequivalent L positions due to the orthorhom- 
bic distortion and/or staggered orbital order, the vari- 
ous mechanisms of polarization are manifested as Fourier 
harmonics in X-ray or neutron scattering probes. Assum- 
ing a general conical structure m r = m,Q + mi cos Q r 
m 2 sin Q r, we obtain the nontrivial Fourier components 
p (l) = Er.e P r + e/2^ M '' (r+e/2) shown in Table I, with 
m± = mi ± im2 and P r + e /2 given in Eq. (1). To com- 
pensate the electrostatic energy loss, each L ion tends 
to shift towards the induced dipole moment and cause 
the lattice modulation with the same wave vectors. P ms 
gives a collinear lattice modulation along the cluster di- 
rection with wave vectors ire + Q and ire + 2Q, where 
the shift of ire is due to the present of two equivalent L 
sites. It becomes uniform for the special case Q-e = ir/2, 
or period four for magnetism. This case is relevant for 
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TABLE I: Wave vectors and directions of the lattice modu- 
lation originating from the mechanisms P ms , P° lh , and P sp . 
P nm always gives a modulation with wave vector ne. 

HoMn03 as pointed out by Sergienko et al. recently [14]. 
Under the spiral (conical) magnetic structure, P orb yields 
a spiral, even noncoplanar, modulation, while P sp gives a 
uniform (conical) modulation. p orb is found only in the 
unfulfilled t2 g system [11]. These results can be used to 
classify the mechanisms of the magnetoelectric couplings 
using X-ray and neutron scattering experiments. 

Equation (1) can be derived from the TM-L-TM clus- 
ter model [8, 11] composed of the L p orbitals and the two 
TM d orbitals which are hybridized with L through the 
Slater-Koster parameters V p da and V p d-w- We assume an 
effective Zeeman field (U/2)m a (a = l/r for the left/right 
TM ion) with an energy separation U , which originates 
from the local Coulomb repulsion and the Hund coupling 
in the magnetically ordered phase. The local energy lev- 
els of the TM d a and L p orbitals are Ed a and E p , respec- 
tively, while the spin-orbit couplings for d and p orbitals 
are denoted A^ and X p . 

The mechanisms responsible for the spin-induced 
dipole moments through spin-orbit coupling can be clas- 
sified as in Table II. Spin-orbit interaction X t2g within 
the degenerate t2 g manifold can be a source of polariza- 
tion, P sp and P orb [8, 11], when the t 2g orbital degrees 
of freedom are unquenched, e.g. V^ g with n = 1,2,4,5. 
It was derived in a previous paper [11] that P olh and 
P sp appear in proportion to min(A t j/V^ £ i 7r , I^V^tt/A) 
and (Ad/A)(V^j7r/A) 3 , respectively, in the t\ g configura- 
tion, where A = E t2g + U/2 — E p . The t% g configuration 
is related to t\ g with the change in the charge transfer 
energy to A = E t2g — U/2 — E p , and the t\ g to t\ g in the 
same manner. Another symmetry operation can be used 
to show that P sp and P orb for t\ g is related to those of t\ g 
by a sign change, as well as those for t\ to t\ g . Due to 
the strong Hund coupling, the t\ g and t\ g configurations 
represent the orbital quenched case with no polarization 
arising from orbital deformation. 

The spin-orbit interaction within the e g manifold van- 
ishes and has no impact on spin-polarization coupling. In 
this case, the spin-orbit interaction A p at the ligand oxy- 
gen site can mediate the spin-polarization coupling P sp , 
as we demonstrate below. In addition, mixing of an occu- 
pied e g level and an unoccupied t2 g , or of an unoccupied 
e g and an occupied t2 g , can occur due to the spin-orbit 
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TABLE II: Classification of the spin-orbit interactions and the 
corresponding types of induced polarization and the relevant d 
electron configuration. Nonzero net magnetic moment for the 
d electrons, S 0, is required to generate magnetic ground 
states. A t2s _ es (A p ) effects decrease (increase) as the electron 
number n increases. 



interaction A t2B _e B in the given TM ion and contribute to 
the polarization P sp . Such mechanism was demonstrated 
for TbMn03 recently [12]. Both mechanisms described 
in this paragraph should apply for arbitrary d n config- 
uration, as long as the net magnetic moment remains 
nonzero to guarantee magnetic ordering. 

A low-n transition metal involving Ti or V will be the 
likely candidate to observe i2 9 -mediated P sp and P orb . A 
spin-canted phase in YVO3 was found recently [15]. For 
such a state we argue that P orb , which is much larger 
than P sp on the atomic scale, is detectable in the micro- 
scopic X-ray or neutron probes. The specific dependence 
of P orb on the local magnetic structure given in Eq. (1) 
and in Table I can be checked experimentally. 

Recently, a number of e g systems with filled t 2g came 
to be classified as a multiferroic: a nickelate N13V2O8 [7] 
with d 8 , and a cuprate LiCuVCU [18] with d 9 . We 
propose that the oxygen spin-orbit-mediatcd mechanism 
should play an important role for these materials. For 
instance, the d 9 Cu sites can be modeled as d x 2_ y 2 01- 
bitals at the TM sites in the cluster model. Taking 
-Vi = V r = V pda [W], P sp is given by[17] 



f« = -2v^„ + L J „)^(^ 



(2) 



where A = E Bg +U/2 — E p , and E Cg is the energy of 
the relevant e g orbital. Here L a , z denotes the integral 
L a . z — (d a \z\p z ) with a — r,l. The d 8 contains two e g 
orbitals with S = 1, that can be decomposed as d 3x 2_ r 2 
and d y 2_ z 2 in the degenerate case. Of these, the latter 
does not couple to the oxygen orbitals and one is left 
with a single d 3x 2_ r 2 orbital per site with the polarization 
again given by Eq. (2). No p orb term is found for pure 
e g models with oxygen spin-orbit interaction. 

The d 7 case with a single e g electron is subject to the 
Jahn- Teller distortion and the lifting of the orbital degen- 
eracy. When the same e g orbitals are occupied for the 
adjacent TM sites, one again has Eq. (2). If not, the ap- 
propriate cluster model should involve different orbitals 
for the right and left TM's, and the inversion symmetry 
breaking yields non-zero values of 



p- = (L r;X -L, ;X )^, P - = I( i;iX _i riX )(^) 3 , 

(3) 

with L ajX = (d a \x\p x ). On the local scale both terms, 
being independent of spin-orbit interaction, are greater 
than P^ p . Except for special circumstances, however, 
pnm an( j pms are nn j^ c _g modulations with zero macro- 
scopic average. Even for the d 9 and d 8 cases, the or- 
thorhombic distortion displaces the oxygen position away 
from the middle of the TM ions, then one has the appear- 
ance of Eq. (3) in the polarization. 



As a concrete example of the A t 



mechanism 



at work, we discuss the case of $i g e g configuration 
found in TbMn03 with Mn 3+ ion. Because of the 
large Jahn- Teller effect, TbMn03 exhibits the rod-type 
d 3x 2_ r 2/d 3y 2_ r 2 orbital ordering shown in Fig. 1(b) far 
above the room temperature[19]. Accordingly we assume 
that the d 3x 2_ r 2 and d 3y 2_ r 2 orbitals are occupied at 
the I and r sites of the TM-L-TM cluster. Assuming 
Vi = V r = —Vpda, one obtains polarization of the spin- 
current form [17] 
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where A c / is the crystal field gap E Cg — E t2 , Pjt is 
the Jahn- Teller energy splitting, and L' z = (di iZX \z\p x ) = 
(d r ,zx\z\p x )[i2, 17]. A similar term is obtained for uni- 
form orbital order, i.e. same d-orbitals for both TM ions 
in the cluster. No orbital polarization P orb is induced 
from the mixing [17]. 

In the rest of the paper, we compare the experimen- 
tal situation with the theoretical analysis for TbMn03, 
which probably presents the most complex situation with 
orbital ordering, unfulfilled t 2g and e g levels, and t 2g —e g 
mixing. The two terms, Eq. (4) and Eq. (2), combine 
constructively to produce P sp , the uniform polarization. 
Quantitative agreement with the measured uniform po- 
larization value in TbMn03 is obtained as follows. 

Matrix elements needed to evaluate Eq. (2) and Eq. 
(4) are calculated by taking the TM-0 bond length 
1.9A and the Clementi-Raimondi effective charges [21] 

Z Mn3d = 10 - 53 and Z 02 P = 4 - 45 - 0ther parameters 
are chosen as Vi = 2V r = —V p d a = —1.2 cV, U = 3 
eV and E Cg — E p = 2 eV. The spin-orbit coupling for 
the oxygen 2p and TM 3d orbitals are chosen from 
X L = 2AE L /(2L + 1), where L = 1,2 for p and d or- 
bitals. The spin-orbit energy splitting is given by AE p = 
37 mcV (A p = 25 mcV), and AE d = 120 mcV (X d = 48 
meV), respectively [20]. With these values we obtain the 
uniform polarization along the c axis P^ p ~ 130 /xC/m 2 
from the oxygen spin-orbit interaction, and Pt2 g -e g ~ 860 
/iC/m 2 from the t2 g — e g mixing with A c f — 2 eV and 
Ej T = 1 cV. The net value P sp - 990 fiC/m 2 multiplied 
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with \m r x m r+e \ w sin(0.287r) in Eq. (1) gives the uni- 
form c-axis polarization of ~ 760 /xC/m 2 , in reasonable 
agreement with the experimental value in TbMn03[4, 5] 
(^700 /j,C/m 2 at 10K). The direction of the polarization 
is also consistent with the experiment. 

The X-ray intensity of the 2Q lattice modulation in 
TbMn03 increases monotonically between the Neel tem- 
perature Tjv for the collinear magnetic phase, and Tl, the 
latter signaling the onset of ferroelectricity, as the tem- 
perature is lowered. Then the signal decreases with lower 
temperature in the non-collincar, ferrroclcctric phase, 
T < Tl. Both features can be understood by taking 
the view that the non-uniform component arises from 
the magnetostriction term, due to the breaking of the 
inversion symmetry in each TM-L-TM cluster from the 
orbital ordering and the orthorhombic distortion. 

The spin configuration in the magnetic phase T < Xjv 
can be described by the ordered spin moment, (sj) = 
S(b cos 9i + c\fl — e 2 sinflj). It evolves from collinear 
(e = 1) to nearly circular (e = 0) pattern with de- 
creasing temperature [6]. In the collinear phase, one has 
e = 1 while the magnitude S increases monotonically 
at the lower temperature. The TM-L-TM cluster cal- 
culation can be generalized to handle these cases, both 
collinear and elliptic, by writing a self-consistent equa- 

tion (sj) = with a weight w. k — exp[/3j(sj) -m,-J. 



Average over the 0(3) vector determines the appro- 
priate Pi. The dipole moment in the collinear/elliptic 
spin configuration can be obtained by taking the average 
of Eq. (1) over the spins m r and m r+e with weights w r 
and uv+ e , for the magnetic period ~ 0.287T configuration 
of the low-temperature TbMnC>3. 

As shown in Fig. 2 (d), in the collinear phase, (e = 1), 
the 2Q intensity arising from the magnetostriction P ms 
monotonously increases with the ordered spin amplitude 
S. Then in the non-collinear phase (e < 1), the 2Q 
intensity arising from P ms decreases for smaller e cor- 
responding to a lower temperature (Fig. 2 (c)). The 
direction of the lattice modulation in TbMnC>3 has not 
been clearly determined yet [22]. A proper identification 
of the lattice modulation direction in the future can be 
used to discriminate the different scenarios of magneto- 
electric coupling. 
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FIG. 2: (color online) (a) Elliptic spiral spins (red arrows) 
with the ellipticity e = 0.6 and the associated local dipole mo- 
ments arising from oxygen spin-orbit interaction are shown as 
P| p (blue arrows, enlarged by 10 3 times) and P™ s (black ar- 
rows) for the observed spin modulation period (Q » 0.287r) as 
in TbMnOa at low temperature. Two types of Mn sites exists 
(blue and white) due to orbital ordering. Statistical weights 
are introduced to account for elliptic spins. P™ s alternates its 
direction due to staggered orbital order, (b) Induced dipole 
moments P^ and |P™f| at the oxygen site Oi in TMi-Ov 
TMi+i cluster in units of 10 _3 L Z and 10~ 3 L X , respectively, 
and Ps( x ) is defined by Li tZ t x \ . (c,d) Main Fourier components 
of the dipole moments, at q — 2Q for |P™j S | and at q = for 
P^, as functions of e for T < Tl in (c) and S with e = 1 
for T L < T < T N in (d). Note that e in (c) and 1/2 - S in 
(d) decrease with lowering temperature. This figure is worked 
out for the orbital ordered model of TbMnOs (Figs. 1) with 
(U, A cf , E JT , Vpda , X P , A d ) = (3.0, 2.0, 1, 1.2, 0.025, 0.048) eV. 
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